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Abstract A simple analytical model for a turbulent flow is proposed, which considers the flow as a 
collection of localized spatial structures that are composed of elementary "cells" in which the state 
of the particles (atoms or molecules) is uncertain. The Reynolds number is associated with the ratio 
between the total phase volume for the system and that for the elementary cell. Calculating the 
statistical weights of the collections of the localized structures, it is shown that as the Reynolds number 
increases, the elementary cells group into the localized structures, which successfully explains the onset 
of turbulence and some other characteristic properties of turbulent flows. It is also shown that the basic 
assumptions underlying the model are involved in the derivation of the Navier-Stokes equation, which 
suggests that the driving force for the turbulent transition described with the hydrodynamic equations 
is essentially the same as in the present model, i.e. the tendency of the system to occupy a statistically 
dominant state plays a key role. The instability of the flow can then be a mechanism to initiate the 
structural rearrangement of the flow to find this state. 
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1 Introduction 

The transition from laminar to turbulent fluid motion occurring at large Reynolds numbers [1,2.3,4, 
[5] is generally associated with the instability of the laminar flow, and this viewpoint is well supported 
by the analysis of solutions of hydrodynamic equations [H0G1]. On the other hand, since the turbulent 
flow characteristically appears in the form of eddies filling the flow field, the tendency to occupy 
such a structured state of the flow cannot be ruled out as a driving force for turbulent transition. 
To examine this possibility, we propose a simple analytical model for the flow and show that as the 
Reynolds number increases, the state of the flow in the form of collection of localized spatial structures 
(eddies) becomes statistically more favorable than the unstructured state, which successfully explains 
the onset of turbulence and some other general properties of turbulent flows. We also show that the 
basic assumptions underlying the model are involved in the derivation of the Navier-Stokes equation, 
which suggests that the driving force for the turbulent transition is basically the same as in the present 
model, i.e. the tendency of the system to occupy a statistically dominant state plays a key role. The 
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instability of the flow at high Reynolds numbers can then be a mechanism to initiate the structural 
rearrangement of the flow to hnd this state. 



2 Model 

Let a fluid flow be represented by a system of N identical particles (atoms or molecules). Assume that 
the particles can form localized structures of Ni particles (with the number of such structures being 
Mi), which are similar to turbulent eddies in that the particles execute an overall concerted motion. 
Since the concerted motion should break down at the microscale level, each i structure can be divided 
into Ni/n elementary "cells" (each of n particles), in which the states of the particles (positions and 
velocities) are uncorrelated. Then, considering the interaction between the localized structures to be 
small and taking into account that the permutations of the particles in the elementary cells, those 
of the elementary cells in the localized structures, and the permutations of the localized structures 
themselves do not lead to new states, the total number of states can be estimated as for an ideal-gas 
system [7] 

m 



Hi [(n!)^/«(7V i /n)!] M *M i ! 



or, with the Stirling approximation x\ w (xje) x to be applicable 

r= (N/n) N e N ^ 

n i (A r l /n) Ar ' M -/"(M l /e)^ 1 ' 

where £\ A^M, = N. 

With Eq. ([lj, it is possible to calculate the most probable size distribution of the structures, 
Mi = Mi(Ni), which maximizes the number of states r and, correspondingly, the entropy S = hxT. 
This distribution reduces the phase space of the system to a subspace of much lower dimension, similar 
to how, according to the dynamical system theory, the phase space shrinks when the system reaches the 
(strange) attractor in the transition to turbulence 4,8. Applying the Lagrange multiplier method to the 
entropy functional, i.e. varying the expression ^\ [—(MiNi/n) \n(Ni/n) — Mi ln(Mi/e)]+a' MiN+ 
j3' J^i MiEi with respect to Mj at two conservation conditions J^- MjA/j = N and MiEi = E [Ei 
is the kinetic energy of fluctuations in i structure, and E is the total kinetic energy), one obtains 

^ = " e Q ' w *+^ (2) 

where and a' and /3' are the Lagrange multipliers. 

The dependence of J3j on Ni is specific rather than universal; it is clearly different, e.g., for homo- 
geneous and shear-layer turbulence. To be definite, assume that the velocity of fluctuations increases 
with the distance r as in the Kolmogorov theory of turbulence for the inertial range of scales [41I51ITU] , 
i.e. 

v(r) - r h (3) 
where h = 1/3. Then, the kinetic energy of fluctuations per unit mass is e(r) ~ r 2h , and ~ 
Jo' e { r ) r2 dr ~ r 2h+3 ~ j\[ 2h / 3+1 _ Also, it is convenient to pass from N to the number of the elementary 
cells in the localized structure qi = Ni/n. With these changes, Eq. ^ becomes 

Mi = q7^e a ^ h/3+1 (4) 

where a and /3 are new constants that are determined from the equations of conservation of the total 
number of particles and kinetic energy 

N/n = J2 M i9i (5) 

i 

and 



E/n 2h / 3+1 =Y, M ^ h/3+1 ( 6 ) 
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Given the values of n, N and E, the constants a and j3 can be calculated by generating a set of qi 
and solving Eqs. ^ and (El with Mj substituted from Eq. Another possibility is to vary a and /3 
to obtain N/n and E/n 2h ' 3+1 as functions of a and /?. 

The parameter N/n can be associated with the Reynolds number Rez, = VL/u, where V and L are, 
respectively, the velocity and linear scale characterizing the system as a whole, and v is the kinematic 
viscosity of the fluid. Since the particles are identical, the 6iV-dimensional phase space is reduced to the 
6-dimensional single-particle space. Then, the linear size of the elementary cell, in which the behavior 
of the particles is uncorrelated, can be taken as the product of the characteristic velocity and length at 
which the state of the particles is uncertain. For a gas fluid, it is the product of the molecular thermal 
velocity c and the mean free path A [TT]. Then, the phase volume of the cell is (cA) 3 , with the number 
of particles in the cell being proportional to this volume, n ~ (cA) 3 . Correspondingly, the total number 
of particles N can be considered to be proportional to the total phase volume for the system (VL) 3 , 
i.e. N ~ (VL) 3 (the density of the fluid is assumed to be constant). Taking into account the kinetic 
theory expression for the gas viscosity v ~ cA (e.g., Ref. [12])) one obtains 

N/n ~ (VL/vf = Re| (7) 

A similar relation is valid for a liquid fluid. Here, the linear size of the elementary cell can be taken 
to be c 2 t, where c 2 presents the fluctuations of the (molecular) kinetic energy per unit mass, and t is 
the mean residence time of the molecule in a settled state. Since the viscosity of liquid is v ~ c 2 r (e.g., 
Ref. [13]), we arrive at Eq. (7k again. In what follows, we will assume for simplicity that the coefficient 
of proportionality in Eq. |7| is equal to 1, i.e. 

Re L = (N/n) 1 ' 3 (8) 



3 Results and Discussion 

3.1 Turbulent Transition Region 

Figure [T] shows a characteristic dependence of the average number of elementary cells in the localized 
structures (q) = Qi^i/ X)j % on the Reynolds number Rez,; the number of elementary cells 
varies from 1 to q max . Points labelled with crosses represent all combinations of a and f3, with these 
parameters varied independently within a m ; n < a < a max and /3 m - ln < < /3 max . However, not every 
combination leads to a physically reasonable bell-shaped distribution [T3] (see also Fig. [i] below) . In 
small fraction of the combinations (typically at large values of /3), Mi does not vanish at q — > g max . 
The fraction of such "wrong" combinations is generally within several percentages (in the given case, 
below 2 percent). Triangles show the dependence (q) on Rez, with those "wrong" a/fi combinations 
excluded. The increase of g max does not change the distribution, except that the range of variation of 
(q) and Rez, is extended to larger values of these quantities. 

It is seen that as Rez, exceeds some characteristic value (Rez ~ 10 2 in Fig.[l]), (q) rapidly increases, 
i.e. the elementary cells group into localized structures, signaling that the flow becomes turbulent at 
these Reynolds numbers. It is noteworthy that the dependence of (q) on Rez, is not unique, i.e. the 
same values of (q) are observed in a broad range of Rez . This indicates that the turbulent state is not 
solely determined by the Reynolds number; rather, as is well-known, it is flow specific, depending on 
the type of the flow, the inlet conditions, the flow environment, etc. [3l|4j[5l[6j. The present statistical 
model is too simple to take these effects into account, but it offers an estimate for the range of Re^ at 
which the turbulent motion can be expected, i.e. the lower and upper bounds of the turbulent transition 
region. The lower and upper boundaries of the manifold of the "correct" points in Fig. [l] (labelled with 
triangles) are determined, respectively, by the maximum values of j3 and a at which the correct points 
are obtained (Fig. [2]). As these parameters increase, the lower boundary shifts to smaller values of 
Re^, and the upper boundary to larger values of ReL (Supplementary Material). More specifically, at 
the lower boundary Re^ ~ exp(— 0.15/3(g)), and at the upper boundary Re^ ~ exp(0.06a(q)). Figure 
[3] shows these dependencies for (q) = 40, which corresponds to approximately the maxima of the size 
distributions of the localized structures (see Fig. [4 and its discussion below). It is remarkable that the 
lower boundary moves until Rez, ~ 10 2 is reached; after that it "freezes" (Supplementary Material). 
The model thus predicts that the flow should be laminar at Rez < Re^ (Re^ ~ 10 2 in the present 
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Fig. 1 (Color online) Average number of the elementary cells in the localized structures (q) as a function of 
Re_L. The number of the elementary cells qi varies from 1 to g max , and parameters a and /3 vary independently 
from a m in = /3 m i n = —1.5 to a m ax = /?max = 1-5 with the interval 0.015. Crosses represent all combinations of a 
and P at g max = 100, and the triangles stand for the combinations which lead to the physically reasonable bell- 
shaped distribution at different g max ; the magenta, blue, red and green triangles are for g m ax = 100,200,300 
and 400, respectively. 
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Fig. 2 The values of the parameters a and /3 at the lower and upper boundaries of the transition region of 
Fig. [1] g max = 200. 



model), and it can remain laminar up to very large values of Re^,. This prediction is in good general 
agreement with the experimental and simulation results [21HJGI]. I R particular, for the pipe flow, in 
his seminal work Reynolds [T] has found that the critical Re number varied from Re^ « 2 x 10 3 to 
Re^ « 1.3 x 10 4 depending on the inlet conditions. The lower bound of the stability of the laminar flow, 
which is observed at large disturbances of the inlet flow, has been confirmed in many experimental 
and theoretical works (e.g., Refs. [15] and [16], respectively), and the upper bound, which is achieved 
in carefully controlled conditions, has been found as high as Re^ ~ 1 x 10 5 [TT]- In the present model, 
the above mentioned lower bound of the stability (Re/, ~ 10 3 ) corresponds to /3 max « 1.25, and the 
upper bound is unlimited (Fig. 3). 



3.2 Other Characteristic Properties of Turbulent Flows 

The model is also consistent with some other characteristic properties of turbulent flows: 

• The dissipation scale 77 obeys the known equation in the Kolmogorov theory of turbulence 



r,/L ~ Re" 3/4 



(9) 
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Fig. 3 The Reynolds number Rez, corresponding to the lower (open triangles) and the upper (solid triangles) 
boundaries of the transition range as a function of parameters /3 and a, respectively; (q) = 40. The dashed 
lines are the best exponential fits, see the text. 



Indeed, similarly to (N/n) 1 / 3 in Eq. M, the quantity q 1 / 3 — (Ni/n) 1 / 3 can be considered to be the 

— 1 /3 

Reynolds number of i structure, i.e. Rei = q i . According to the Kolmogorov theory Re^ ~ 1 
characterizes the dissipation scale n. Let us rewrite qi — Ni/n as (Ni/N) x (N/n) and take into 
account that v t /V ~ [h/L) 1 / 3 [Eq. so that Ni/N ~ [vik/{VL)} 3 ~ (h/L) A . Then, 

q t . ~ {N/n)(k/L) 4 = Rei(Z,/L) 4 (10) 

and the equality qi ~ 1 (with ij = 77) leads to Eq. (p}]). This, in particular, suggests that the linear 
size of the elementary cell could also be determined as rjv v , where v v is the velocity increment at the 
dissipation scale. 

• The probability density distributions of the linear sizes of the structures Mi(li/rj), FigjiJ are 
in general agreement with the results of the direct numerical simulation of isotropic homogeneous 
turbulence by Jimenez et al. [14] . who found that the vortex radius distributions are bell-shaped, 
practically do not shift with Re r, (with the radius r measured in rj units), and have maxima at r ss 3rj. 
According to Eqs. ^ and (10), k/rj rj (qi) 1 / 4 , so that the maxima of the distributions at k/rj ~ 2.5 



correspond to q t w 40. 

• The structure functions and energy spectrum are similar to those in the Kolmogorov-Obukhov 
theory [91IT8]. Since the dependence v = v(r) is assumed to be the same for all localized structures, the 
structure functions are independent of the size distributions of localized structures. With Eq. the 
structure function of p-th order is S p (f) = (Svf) ~ l p / 3 , where I is the space increment, and the energy 
spectrum is E(k) ~ fc~ 5 / 3 , where k — 2tt/1 is the wave number. 



3.3 Connection with the Hydrodynamic Equations 

The above consideration shows that the proposed model captures characteristic properties of turbulent 
flows, in particular, those obtained with the hydrodynamic equations. Then, the natural question to 
ask is: Why do two apparently different approaches, the hydrodynamic equations and the present 
model, lead to similar results? Or, what is common to these approaches to lead to similar results? 
The present model is based on three assumptions: i) the particles are identical, ii) there exists an 
elementary cell of the phase space in which the state of the particles is uncertain, and iii) the kinetic 
energy of the fluctuations is as in the Kolmogorov theory for the inertial interval of scales. The latter 
assumption seems to be not very essential for our purpose. For example, even if a quasi-solid rotation 
of the localized structures (eddies) is supposed, which highly overestimates their kinetic energy because 
the absence of the decrease of the velocity outside the eddy core Q]5], the statistical preference of the 
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Fig. 4 Characteristic linear size distributions of the localized structures, 1 < Qi < 200. Solid lines show the 
Mi(li/rf) distributions averaged over all combinations of a and /3 for which ReL are close to those indicated 
at the curves within 1 percent (typically one of such distributions dominated). The dashed curves are the 
Gaussian fits to the distributions. 

structured flow is still preserved (Supplementary Material). The attention should thus be focused on 
the former two assumptions. As well known, the Navier-Stokes equation for a gas fluid can be derived 
from the Boltzmann equation (the Chapman-Enskog expansion |12|). which, in turn, can be obtained 
from the Liouville equation (see, e.g., Ref. [M]). In the transition from the Liouville equation to the 
Boltzmann equation, first, the many-particle distribution function is reduced to the single-particle one 
under the assumption that the particles are identical (the Bogoliubov-Born-Green-Kirkwood-Yvon 
hierarchy, see, e.g. Ref. [H]), and, secondly, because the gas density is taken to be low, the interaction 
between the particles is assumed to be a random two-body collision, so that the positions and velocities 
of the particles are uncertain within the mean free path A and the thermal velocity c, respectively. It 
follows that both the assumptions are implicitly present in the Navier-Stokes equation [22) . Therefore, 
although the hydrodynamic (Navier-Stokes) equations are incomparably rich in the description of the 
turbulent phenomena, because they are able to take into account the specific conditions under which 
the transition takes place (the inlet conditions, flow environment, etc.) and give a dynamic picture of 
the transition, the driving force for the transition can have the same statistical origin as in the present 
model. 



4 Conclusions 

A model for turbulent fluid flow has been proposed that considers the flow as a collection of local- 
ized spatial structures (eddies) and estimates the statistical weights of these collections for different 
Reynolds numbers. The Reynolds number is associated with the ratio between the total phase volume 
for the system and that for an elementary cell in which the state of the particles is uncertain. It has 
been found that i) the model successfully explains the onset of turbulence as well as some other charac- 
teristic properties of turbulent flows, and ii) the basic assumptions underlying the model, i.e. that the 
particles are identical and an elementary cell of the phase space exists where the state of the particles 
is uncertain, are implicitly present in the Navier-Stokes equation. Considered together, these findings 
suggest that with all the variety of dynamic phenomena that the hydrodynamic (Navier-Stokes) equa- 
tions describe [5], the driving force for the turbulent transition is basically the same as in the present 
model, i.e. the system tends to occupy a spatial state that is statistically dominant at a given Reynolds 
number. The instability of the flow can then be a mechanism to initiate the structural rearrangement 
of the flow to find this state. 

The present statistical aspect of the turbulent transition is believed to be important for a more 
comprehensive understanding of the nature of this challenging phenomenon as well as of the other 
related problems in the field of pattern formation [23 . 

Acknowledgements I thank B. Ilyushin and D. Sikovsky for useful discussions. 
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Supplementary Material: 

Variation of the transition region with parameters a and j3. Figure [ST] illustrates how the 
distribution of Fig. 1 at q max = 200 changes with variation of a and j3. The lower boundary shifts to 
the smaller values of Re/, , and the upper boundary to higher values of Re^ . 



60 h 




Re L 

Fig. SI Average number of the elementary cells in the localized structures (q) as a function of Re^; 1 < q t < 
200. The cyan triangles are for —1.3 < a < 1.3 and —1.3 < /3 < 1.3, the blue triangles are for —1.5 < a < 1.5 
and -1.5 < /? < 1.5 (as in Fig. 1), and the red triangles are for -2.0 < a < 2.0 and -2.0 < (3 < 2.0. 



Variation of the bounds of the flow stability with the Reynolds number. Figure S2 shows 
how the Reynolds numbers at the lower (a) and upper (b) boundaries of the transition range change 
with a and /3. While the upper boundary monotonically shifts to larger values of Re^ as a increases, 
the lower boundary shifts to smaller values of Re^ until j3 = f3* is reached (/3* w 1.725 at (q) = 20, 
and P* R3 1.5 at (q) = 40). At larger values of f3 the lower boundary "freezes", as, e.g., it is seen from 
Fig. [ST] 
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Fig. S2 The Reynolds number Ret corresponding to the lower (panel a) and the upper (panel b) boundaries 
of the transition range as a function of parameters f3 and a, respectively. The empty and solid triangles are for 
(q) = 20 and (q) = 40, respectively. The dashed lines show the best exponential fits to the data. 



Quasi-solid rotation of the localized structures. If a quasi-solid rotation of the localized 
structures is supposed, in which case Ei ~ N^ 3 , (q) rapidly increases with Re^ as previously, i.e. the 
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Fig. S3 Average number of the elementary cells in the localized structures (q) as a function of Rez, for quasi- 
solid rotation of the structures. 1 < qi < 100, and —0.25 < a < 0.25 and —0.25 < f3 < 0.25. Crosses represent 
all combinations of a and /3, and magenta triangles stand for the combinations which lead to the physically 
reasonable bell-shaped distribution. 



statistical preference of the structured flow is preserved (Fig. S3). The difference is that the transition 
region is located within unreasonably low Re numbers. This is because the kinetic energy of the eddies 
is highly overestimated due to the absence of the decrease of the velocity outside the eddy core, which 
is characteristic of more realistic models, e.g. the Rankine and Lamb-Oseen vortices pQ. 
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